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RBC and UKQCD Collaborations 

We report isovector form factors and low moments of structure functions of nucleon in numer- 
ical lattice quantum chromodynamics (QCD) from the on-going calculations by the RIKEN- 
BNL-Columbia (RBC) and UKQCD Collaborations with (2h-1) dynamical flavors of domain-wall 
fermion (DWF) quarks. We calculate the matrix elements with four light quark masses, corre- 
sponding to pion mass values of m^ ~ 330-670 MeV, while the dynamical strange mass is fixed 
at a value close to physical, on (2.7 fm)^ spatial volume. We found that our axial charge, g^, at 
the lightest mass exhibits a large deviation from the heavier mass results. This deviation seems 
to be a finite-size effect as the gA value scales with a single parameter, nij^L, the product of pion 
mass and linear spatial lattice size. The scaling is also seen in earlier 2-flavor dynamical DWF 
and Wilson quark calculations. Without this lightest point, the three heavier mass results show 
only very mild mass dependence and linearly extrapolate to = 1 . 16(6). We determined the four 
form factors, the vector (Dirac), induced tensor (PauU), axial vector and induced pseudoscalar, at 
a few finite momentum transfer values as well. At the physical pion mass the form-factors root 
mean square radii determined from the momentum-transfer dependence are 20-30% smaller than 
the corresonding experiments. The ratio of the isovector quark momentum to helicity fractions, 
{x)u-d/ {x)Au-Ad is in agreement with experiment without much mass dependence including the 
lightest point. We obtain an estimate, 0.81(2), by a constant fit. Although the individual mo- 
mentum and helicity fractions are yet to be renormalized, they show encouraging trend toward 
experiment. 
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1. Introduction 

We report isovector form factors and low moments of structure functions of nucleon in numer- 
ical lattice quantum chromodynamics (QCD) from the on-going calculations by the RIKEN-BNL- 
Columbia (RBC) and UKQCD Collaborations [|l]-^ with degenerate up and down and a heavier 
strange flavors represented by domain- wall fermions (DWF) [Q-^- 

The isovector form factors are defined in the following two equations: 

{p\VniO)\p) = Up [YnFi{q^) + 0^yqyF2{q^)/2mN] Up, (1.1) 
{p\A^{0)\p) = Up [Y^Y5GA{q^) + iqixY5Gp{q^)\ Up, (1.2) 

where = uj^u — dy^d and = uYnYsu — dy^jsd are isovector vector and axial vector currents, 
respectively. These form factors are experimentally measured in neutron decays and other elec- 
troweak transitions of nucleon All the form factors can be calculated numerically on the lattice 

S|]- 

The structure functions are measured in deep inelastic lepton scatterings off nucleon [^. For 



their definitions we refer the readers to an earlier RBC publication [10| and references cited there 
in. In this report we discuss some of their isovector low moments such as the momentum fraction 
{x)u-d, helicity fraction {x)/^u-M, transversity (l)^^-^^ and twist-3 d\. 

2. Formulation 

The matrix elements corresponding to linear combinations of the form factors eqs.( |l.l[ ) and 



(1.2) are determined from the ratio of the three-point and two-point functions 1 1 1] 

G^{t'-t,q)G^{t-to,0)G^{t'-to,0 



(^'^^ = G,-(.'-.o,0) 
where G^ ' {t,q) is three-point function, 



G^{t' - t,0)G^{t - to, q)G^{t' - to,q) 



(2.1) 



G^,'^{t,q) = \Tr[^{0\x{t',OMt,qmto-q))] , (2.2) 

with the current = V^,A^, the projector ^ = ^-Y^i^t),^^757z{^5z)^ the spatial momentum 
transfer q, and x being the nucleon field, to and t' are the sources of the three-point function, 
and G2 is the two-point function with the point(L) or the gauge invariant Gaussian smearing(G) 
operator sink. All the sources for the two- and three-point functions are calculated with the gauge 
invariant Gaussian smearing operator. R^" [t) at fo <C f ^ f will be constant, and corresponds 
to the linear combination of the form factors multiplied by kinematic factors. The form factor is 
obtained by solving the linear equation. The details of the equations are in Ref. 1^]. 

The axial charge is calculated by the ratio of the three-point function G^' (t) / gJ"^' (t) = 
ZySa/Za at zero momentum transfer. This ratio gives the renomalized axial charge gA, because 
the axial vector and vector currents share a common renormalization, Za = Zy, up to second order 
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Figure 1: An Edinburgh plot obtained from the present and related RBC and UKQCD joint ensembles with 
(2+l)-flavor dynamical DWF quarks. 



discretization error thanks to the well-preserved chiral symmetry of the domain wall fermions. We 
confirm that these renormalization constants agree within 0.5% accuracy in the chiral limit. 

The matrix elements related to the structure functions, such as the momentum fraction, helicity 
distribution, moment of transversity and di , are calculated by the ratio of the three-point function 
to the two-point function, G^'"^ {t)/G2it'), at zero momentum transfer. Definitions of the oper- 
ators ^' are listed in detail in Ref. [|T0|]. At the lightest quark mass rrif = 0.005 the three-point 



function with the temporal direction of the conserved vector current [12| ^ at q = is used for the 
denominator of the ratio instead of G^(f'). This ratio gives same matrix element because a relation, 
Gf'^'it)/G^{t') = 1, is satisfied for < ? < t'. 

3. Numerical setup 



The calculations are performed on the QCDOC dedicated computers [13] at RIKEN-BNL 
Research Center (RBRC) and University of Edinburgh. Descriptions of the ensembles used are 
found in RBC-UKQCD publications [|]-^]. We use a combination of Iwasaki rectangular gauge 
action [14] with the gauge coupling set at 2.13 and domain-wall fermion [Q-^ quarks with the 
domain- wall height set at 1.8. A 24^ x 64 lattice is used. We fix the dynamical strange mass at 
0.04 and generated four different ensembles with degenerate up and down mass each at 0.03, 0.02, 
0.01 and 0.005 in lattice units. From these we estimate the lattice cut off to be about a^^ of 1.73(3) 
GeV [15] with the Q. baryon mass. The 24^ lattice spatial volume thus corresponds to a physical 
volume of about (2.74(4)fm)-'. The residual quark mass that parameterizes the mixing of the two 
domain walls across the = 16 fifth dimension is estimated to be about 0.0031. The physical 
strange mass is estimated as about 0.035(1) plus the 0.0031 residual mass from the squared mass 
ratio of kaon and . 

For the nucleon matrix element calculation reported here, 106, 98, 356 and 178 configurations 
are used respectively for light quark mass values of 0.03, 0.02, 0.01 and 0.005. We see a reasonable 
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behavior in both pion and nucleon mass in their approach to the chiral limit, as shown in Figure 
[I]. We obtained pion masses of m;i = 0.67, 0.56, 0.42 and 0.33 GeV, and nucleon m/v = 1.56, 1.39, 
1.25 and 1.15 GeV, respectively from the ensembles with the quoted light quark mass values. 

Four measurements are carried out for each configuration to improve statistics, with the source 
set at time slices to = 0, 16,32,48 or to = 8, 19,40,51 except at nif = 0.005. In order to reduce 
computational cost at rrif = 0.005, we employ non-relativistic quark field and a double source 
method where the two sources, either at (0,32) or (16,48), are set for one quark propagator. In 
the three-point function the source and sink operators are separated by 12 time slices to reduce 
excited-state contamination as much as possible. 



4. Form factor results 



4.1 Axial charge 

Figure § shows our result of the axial charge normalized by Zy. The (2.7 fm)^ volume data are 
well determined and the statistical uncertainties are less than 8%. The data are almost independent 
of the pion mass squared except the lightest point. The lightest pion mass data is about 15% smaller 



than the other pion mass data. An earlier 2-flavor calculation by the RBC Collaboration [ ]16[ ] with 
the spatial volume (1.9 fm)^ and 1/a = 1.7 GeV showed a similar trend, but with the downward 
behavior setting in at a heavier pion mass than the current 2+ 1 flavor case. 

We suspect that this pion mass dependence driving the axial charge away from the experiment 
at light quark mass values is caused by the finite lattice volume: In general such finite volume effect 
is expected to grow as we set the quark mass lighter as such lighter quarks fluctuate more. This 
interpretation is not inconsistent with the observed behavior of the (2-1-1)- (present) and 2-flavor 
([[l^) DWF results. Furthermore the finite volume effect is larger on smaller spatial volume when 
the quark mass is same. More quantitatively, we observe in the figure that the 2-flavor result from 
the (1.9 fm)^ volume significantly decreases at nijc = 0.24 GeV^, while the (2-i-l)-flavor results from 
the (2.7 fm)-' volume does not even at m| = 0.17 GeV^. The similar behavior was also observed in 
an earlier small- volume, quenched study f^]. 



a N=2+l(2.7fm) 

■ N=2+l(1.8fm) 

^ N=2(1.9fm) 

^ experiment 



1 



1 



mJGeV ] 

Figure 2: Axial charges with 2+1 flavor and 2 flavor. Dashed line represents linear chiral extrapolation of 
larger volume, 2h-1 flavor data without lightest point. 
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Also shown in the figure is a set of (2+l)-flavor results from a smaller volume [17], (1.8 fm)-', 
with the pion mass, spatial volume and lattice spacing comparable to the 2-flavor calculations. 
However, the result at the lightest pion mass suffers from a large statistical fluctuation and prevents 
us from deciding whether there is a similar finite-size effect here. 

In order to compare the results from the (2+1)- and 2-flavor calculations, we plot the axial 
charge against a dimensionless quantity, nij^L, as presented in the left panel of Figure ||. The 
(2-i-l)-flavor, larger volume results and the 2-flavor ones align well with each other, and suggest a 
monotonic dependence on nijiL: in other words, an mj^L scaling. The (2-i-l)-flavor results from the 
smaller volume are also consistent with this m^cL scaling except for the lightest point that suffers 
from large statistical fluctuation. This large statistical fluctuation itself can be another manifestation 
of a large finite-size effect. Nevertheless we plan to improve the statistics of this lightest point so 
as to test the reliability of this m^rL-scaling interpretation. 

We also plot the axial charge calculated with 2 flavors of dynamical Wilson and improved 
Wilson quarks respectively by LHPC/SESAM |jl8|] and QCDSF [^] collaborations against m^L 
(see the right panel in Figure ^.) These calculations were performed at various different spatial 
volumes, pion masses, and the gauge couplings. Like our (2-1-1)- and 2-flavor DWF results in the 
above, all of these were calculated at unitary points where the sea and valence quark masses are 
equal, K^ea = K'ai- These Wilson quark results also seem to suggest a similar scaling in rrijiL, with 
a downward behavior setting in at m^L around and below 6. 

In the above we discussed that the downward shift away from the experiment of the axial 
charge at lighter quark mass values may well be caused by finite lattice volumes. We found the 
axial charge is well described by a monotonic scaling in the dimensionless parameter, mj^L. This 
mj[L scaling seems common among both DWF and Wilson dynamical quark calculations performed 
at unitary points. Therefore, we seem to have a strong case to suspect a large finite-size effect in 
nucleon electroweak matrix elements in the quark mass range relevant for extrapolating to the 
physical or chiral point. 

Thus we carry out chiral extrapolation of the axial charge using only the larger volume results 
and without the lightest pion mass point. We simply use a linear function of the pion mass squared, 
because there now are only three available data points which do not suggest any non-linear behavior 
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Figure 3: Axial charges with our dynamical domain wall(left panel) and dynamical Wilson fermions(right 
panel) obtained by LHPC/SESAM |Jl|] and QCDSF [|l|] collaborations as a function of mj,L. 
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in the pion mass squared. The extrapolation is presented also in Figure^: we obtain gA = 1.16(6) 
at the physical pion mass rrijc = 0.14 GeV. This extrapolated value is 8.3% smaller than the experi- 
ment. Note this deviation may also be caused by a finite volume effect which may be present even 
at the heavier quark mass values used in the estimation. 

We need more detailed study to clarify this possibly large finite volume effect, not only in 
the axial charge but also in other form factors and structure functions. We plan a larger volume 
calculation at the same cut off in the near future. 

4.2 Isovector Dirac form factor 

The left panel of Figure § shows our result of the isovector Dirac form factor at the four quark 
mass values plotted against the momentum transfer squared, q^. The form factor is renormalized 
by Zy, in other words, normalized by 1 /Fi (0). Traditionally the experimental Dirac form factor is 
considered to be approximated well by a dipole form, 

Fl{q^) = \/{l+q^/Mf)\ (4.1) 

where Mi denotes the dipole mass. This traditional experimental dipole fit is shown in the figure 
also, represented by the dashed line with Mi =0.857(8) GeV [^]. Thus it is convenient to fit the 
present, lattice-calculated Dirac form factor by the dipole form as well. The fit results are presented 
in the figure as solid lines for each quark mass value. With a mild dependence on the quark mass, 
there is tendency for the calculated results to approach the experimental line as the quark mass is 
decreased. 

The Dirac root mean square (rms) radius is determined by the dipole mass by a relation: 

^JJf^ = V^/Mi. Thus the experimental value is 0.794(4) fm. The right panel of Figure ^ shows 
our results for the Dirac rms radius obtained by the dipole fit. The pion mass dependence is again 
mild as it was in the dipole fit discussed in the above. While the result at the lightest pion mass is 
trending toward the experiment with a large statistical error, we cannot exclude that large finite-size 
effect as discussed in the above subsection ^ does not affect the data. We nevertheless decided 
to exclude this point from our chiral extrapolation which is represented by the dashed line in the 
figure: we employ a linear chiral extrapolation due to the mild m| dependence of the heavier data. 
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and obtain y (r^ ) = 0.59(4) fm. The result reproduces about 74% of the experimental value at the 
physical pion mass. As can be seen from the figure, the lightest point is consistent with this fit with 
its large statistical fluctuation. 

4.3 Isovector Pauli form factor 

The isovector Pauli form factor is the induced tensor part of the vector current matrix element. 
We can calculate it at the same time with the Dirac form factor except at zero momentum transfer 
where the kinematics prevent us. Figure § presents the results of our calculation of this form 
factor at each quark mass plotted against the four momentum transfer squared. The form factor is 
renormalized by Zy again. 

The results, unlike the Dirac form factor, suffer from statistical fluctuation. Yet they can be 
fitted by the dipole form, 

F2iq^)=F2{0)/{l+qyMj)\ (4.2) 

The fit results are shown in the figure. There is a large quark mass dependence: the form factor 
decreases as the quark mass is decreased. The three heavier quark mass results seem to approach 
the experiment, shown as the dashed line with the dipole mass M2 = 0.78(2) GeV [^. The lightest 
quark mass result is an exception, however. Note it suffers large fluctuation. 

We determine the anomalous magnetic moment ^2(0) and the Pauli rms radius \J~ir^ = 
y/Vl/M2 from the dipole fit. Figure |^ shows the former in the left and the latter in the right panel 
as functions of the pion mass squared. 

The anomalous magnetic moment shows only a mild mass dependence, and is almost consis- 
tent with the experiment, ^2(0) = ;Up — /x„ — 1 = 3.70589 even at the heaviest point {pLp and /x„ are 
the magnetic moments of proton and neutron respectively.) Though we again cannot exclude the 
possibility of large finite-size effect at the lightest point, the result there is almost consistent with 
the linear fit to the rest which gives an extrapolation to the physical point, ^2(0) = 3.0(5). 

In contrast, the rms radius shows a strong mass dependence approaching the experiment as 
the pion mass decreases. Again we cannot exclude the possibility of large finite-size effect at 
the lightest point. Omitting the lightest point the linear fit gives an extrapolated rms radius of 
Y'' {r\) = 0.69(9) fm at the physical point. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 o.s 
q^[GeV^] 

Figure 5: Isovector Pauli form factor as a function of momentum transfer squared. 
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Figure 6: Pauli form factor at zero momentum transfer and Pauli rms radius are presented in left and right 
panels, respectively. Dashed lines are chiral extrapolations without lightest pion mass data. 



These extrapolated values of anomalous magnetic moment and Pauli rms radius are about 20 % 
lower than the corresponding experiments. While the former almost catches the experiment within 
one standard deviation, the latter is about two standard deviations away. 

4.4 Isovector axial vector form factor 

In this section we focus only on the momentum transfer dependence of the axial vector form 
factor: We normahze the form factor by its value at zero momentum transfer respectively for each 
quark mass. The left panel of Figure ^ shows the results after these normalizations, GA{q^)/GA{0). 
The experimental form factor is again traditionally considered to be fitted well by the dipole form. 



GA{q^)/GA{0) = l/{l+qVMl)\ 



(4.3) 



with the experiments giving Ma = 1.03(2) GeV pOO for the axial vector dipole mass. The experi- 
mental fit is shown by the dashed line in the figure. 

Fits to the calculations with the dipole form are represented by the solid lines and well describe 
the calculations for the heavier three quark mass values. The lightest quark mass results suffer large 
statistical errors and fluctuations. This is because the presented values are normalized by Ga(0), 
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Figure 7: Isovector axial vector form factor and axial vector rms radius are presented in left and right panels, 
respectively. Dashed line in right panel is chiral extrapolation without lightest pion mass data. 
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a quantity directly proportional to the axial charge which itself suffers large statistical errors and 



fluctuations as discussed in the subsection fl.l| . The mass dependence here is even milder than in 
the case of the vector current Dirac form factor. 

The axial charge rms radius is determined from the dipole mass, y'^^'i)^ ~ V^/M^, and is 
0.666(14) fm in the experiment. The calculated axial charge rms radii from the fits are shown in 
the right panel of Figure ^ plotted against the pion mass squared. While the result increases as the 
pion mass decreases, it is about 30% smaller than the experiment. The lightest pion mass data is 
omitted in the following chiral extrapolation, because we cannot rule out a large systematic error 
stemming from the suspected large finite-size effect. However the result would not change, as can 
be seen from the figure, if we included the point. We carry out a linear fit and extrapolation with the 
heavier three mass values and obtain 0.49(4) fm at the physical pion mass. The result reproduces 
73% of the experiment. 



4.5 Isovector induced pseudoscalar form factor 

The induced pseudoscalar form factor, Gp{q^), is obtained as a part of the matrix element of 
the axial vector current. This form factor is expected to have a pion pole, so its momentum-transfer 
dependence should be different from the other form factors. 

The left panel of Figure ^ shows the calculated Gp{q^) renormalized with Zy as plotted against 
the momentum transfer squared at each quark mass. Note that the values at the smallest are 
almost 8 and much larger than other form factors. In addition, while it may be hard to observe due 
to large errors, the values at the smallest increase as the quark mass is decreased except at the 
Ughtest quark mass. This behavior is consistent with pion pole dominance. 

The induced pseudoscalar form factor is related to the axial vector form factor through the 
so-called partially conserved axial vector current (PCAC) relation which is a manifestation of 
spontaneously broken chiral symmetry. In the traditional PCAC current algebra, the celebrated 
generalized Goldberger-Treiman relation, 

Gp{q^) = 2mNGA{q^)/{q^+ml), (4.4) 

is obtained at low q^. The denominator on the right-hand side of this relation corresponds to the 
pion pole. We investigate the relation in our results through a quantity, {q^ +m\)Gp{q-^) / {2M!^GA{q^))- 
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Figure 8: Isovector induced pseudoscalar form factor and ratio of pseudoscalar form factor to axial vector 
form factor are presented in left and right panels, respectively. 
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Figure 9: Left panel is the ;rA^A^ coupling evaluated by GT relation with experimental and measured f^. 
Right panel is induced pseudoscalar coupling for muon capture determined with generalized GT relation 
and induced pseudoscalar form factor Dashed lines are chiral extrapolations without lightest pion mass 
data. 



If the relation holds we obtain unity for this quantity. The right panel of Figure || shows the quan- 
tity in our calculation stays close to unity, and has no significant dependence. We can simply 
fit these results by a constant for each quark mass respectively, whose results are presented in the 



figure as well. All the fit results are consistent with the experiments [21-23] within the larger error 
of the experiments. 

4.6 TiNN coupling 



The original Goldberger-Treiman (GT) relation [24], 



gnNNfn = '^-iriNgA, 



(4.5) 



equates a combination of quantities at the pion pole, the kNN coupling, g^NN, and the pion decay 
constant, /;;:, on the left with another combination of quantities at almost zero momentum transfer, 
the nucleon mass and axial charge. As such it suffers from the mismatch in momentum transfer if 
we substitute the experimental values for the quantities, such as fj^ = 93 MeV, niN = 940 MeV and 
gA = 1.269, to obtain a value for the kNN coupling, g^NN- Thus it becomes interesting how much 
better or worse our lattice calculation does in this regard. 

In the left panel of Figure ^ we show two such calculations each for the nNN coupling, g^NN- 
one uses the experimental value of fj^ and another the lattice-calculated values at each quark mass, 
plotted against the pion mass squared. The results with the experimental fj^ exhibit a significant 
slope in terms of while that with the calculated fji is almost fiat. In both methods the results 
at the lightest mass show significant downward shift away from the trend set by the three heavier 
mass values. This of course is another manifestation of the large deviation observed in the axial 



charge which was discussed in detail in subsection [4.4 For the chiral extrapolation we simply 
employ a linear fit form and exclude the lightest mass point. We obtain the results at the physical 
pion mass, gjiNN = 13.2(9) with the experimental /;;: and g,tNN = 12.1(7) with the calculated fj^. 
The extrapolated results at the physical pion mass are consistent with the experiment obtained from 



forward K-N scattering data [25] 



10 



Nucleon form factors and structure functions with Nf=2+1 dynamical DWF T. Yamazaki and S. Ohta 



4.7 Muon capture 

The induced pseudoscalar coupling, gp = m^Gp{ql), is defined with the muon mass and 
the induced pseudoscalar form factor Gp at the momentum transfer squared, = GeV^. 
This is for convenience in application for muon capture, p + jJ.^ —^n + v^, where the two-body 
nature of the process defines the momentum transfer, q^. 

Gp{q^) can be described by the pion pole dominance form through the generalized Goldberger- 
Treiman relation (^^, and a constant a which corresponds to the difference from unity of the 
quantity, {q^ + m\)Gp{q^) / {2M]siGA{q^)), at each quark mass. These were summarized in Figure 
Thus, gp is determined by 

IniNGAiql) 2mN gA 



where we use the dipole form of the axial vector form factor (fl.3[). In order to subtract the strong 
pion mass dependence stemming from the pion pole, we use the physical pion mass in the pion 
pole. On the other hand we use the calculated values for mj^, gA, Ma and 01. The right panel of 
Figure ^ shows that the result denoted by circle is almost linear as a function of the pion mass 
squared and decreases toward the experiment for the three heavier mass values. Again the lightest 
mass result is an exception: We suspect the significant drop here at the lightest point is caused by 
the finite volume effect in gA as discussed in Sec. p~T| . 

We also determine gp from Gp{q^) directly with the subtraction of the pion pole. At each 
quark mass the quantity Gp{q^){q^ + m\) is fitted by the dipole form, and then extrapolated to <7^. 
In the figure the extrapolated result to q^ normalized by the pion pole with the physical pion mass 
is presented by square symbol. The result from Gp{q^) has larger error, whilst the two results agree 
within the statistical errors at each pion mass. We carry out chiral extrapolations with a simple 
linear function of the pion mass squared without the lightest mass. The results at the physical pion 



mass are consistent with each other, and also with the recent experiment [ ]21[ ] and analysis 

5. Low moments of the structure functions 

5.1 Quark momentum and helicity fractions 

Let us first discuss the naturally renormalized ratio of the isovector quark momentum fraction 
{x)u-d to helicity fraction {x)/^u-Ad- Since the two fractions are related with each other by a chiral 
rotation, they share a common renormalization. And because the DWF quarks preserve the chiral 
symmetry very well, as parameterized by our small residual mass of 0.0031, the ratio is naturally 



renormaUzed on the lattice. Our results are summarized in Figure |10[ The ratio does not show any 
appreciable dependence on the quark mass, albeit with large statistical errors, and is in agreement 
with the experimental value. The constant fit result with all four pion mass data, {x)u-d / {^)au-m = 
0.81(2), is also consistent with the experiment. 

Each of the fractions itself, unlike the ratio, is yet to be renormalized, but suggests interesting 



behavior as shown in Figure 1 1 . They both trend down toward the experimental value at the lightest 
quark mass, after staying almost constant for the three heavier mass values. This may well be 
related to the finite-size effect we suspect for the form factor values, but not necessarily so as the 



11 



Nucleon form factors and structure functions with Nf=2+1 dynamical DWF T. Yamazaki and S. Ohta 



0.9 


1 i i i i 1 

" experiment 




' 1 ' 


i i i 1 










( 


» 


( 


[rrrrrrrJ 


0.8 




i 










0.7 










<x> ^ 


0.6 




... 1 




. 1 . 


Au-Ad J 
. . . 1 . . . . 1 . . . . ■ 



m^[GeV^] 

Figure 10: The naturally renormalized ratio of the isovector quark momentum fraction to helicity 

fraction {x)Au-Ad- 



Structure functions probe different, deep inelastic, physics from the elastic form factors. We will be 
better able to discuss these quantities after finishing the lattice non-perturbative renormalizations 
for them in the near future. 



5.2 Transversity 

The transversity, {'\-)su-Sd^ can be measured by the RHIC Spin experiment in the near future. 
In the present calculation it shows a similar behavior with the quark momentum and helicity frac- 
tions: it stays almost constant for the three heavier quark mass values and then trends down at the 



Ughtest mass (see the left panel of Figure |12[) This quantity is also yet to be renormalized, but will 
soon be, and then will provide a prediction for the experiments planned in the near future 



5.3 Twist-3 di moment 

The di moment of the twist-3 part of the polarized structure function, according to Wandzura 



and Wilczek [27], is small when it is calculated perturbatively. It need not be small under a non- 
perturbative, confining environment within the nucleon. Our calculation for this di moment is 



summarized in the right panel of Figure |12[ The result suggests, though again not renormalized 
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yet to be renormalized. 
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yet, that it is small, consistent with Wandzura and Wilczek. We note also that the results at the 
lightest quark mass show some deviation from the linear extrapolations from the three heavier 
mass values. 

6. Conclusions 

We calculated the isovector form factors and low moments of structure functions of the nucleon 
with Nf = 2+\ dynamical domain wall fermions at 1.7 GeV cutoff on a (2.7 fm)^ spatial volume. 
The axial charge at the lightest quark mass is about 15% smaller than the other mass data, and it 
seems affected by a finite volume effect. It scales with a single parameter, rujiL, the product of pion 
mass and the linear spatial lattice size. We confirmed similar scaling in earlier Nf = 2 dynamical 
DWF and Wilson fermion calculations. Without the lightest point the axial charge is estimated as 
1.16(6) by a linear extrapolation to the physical point. The root mean square radii for the form 
factors except the induced pseudoscalar are determined and at the physical pion mass are 20-30% 
smaller than experiments. The nNN coupling and induced pseudoscalar coupling are found to be 
consistent with experiments. 

We found the renormalized ratio of the isovector quark momentum fraction {x)u-d to helicity 
fraction {x)^^^^^^^ is in agreement with experiment. Their individual values, though yet to be renor- 
malized, show an encouraging trend toward the experimental values at the lightest quark mass. 
Their non-perturbative renormalization will be completed soon. We calculated the bare transver- 
sity which will provide a prediction when its non-perturbative renormalization is completed in the 
near future. We found the twist-3 d\ moment is small, consistent with the Wandzura- Wilczek 
relation. 
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